Spin canting as a result of the competition between stripes and spirals in cuprates 
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Based on the extended Hubbard model we calculate the energy of stripe and spiral ground states. 
We find that uniform spirals get favored by a large t'/t ratio but are unstable at small doping towards 
stripes and checkerboard textures with spin canting. The structure of these inhomogeneities also 
depends on t'/t and the associated spin currents may induce a small lattice distortion associated 
with local dipole moments. We discuss a new kind of stripe which appears as a domain wall of the 
antiferromagnetic (AF) order parameter with a fractional change of the phase of the AF order. For 
large \t [t\ spirals can be stabilized under certain conditions in the overdoped regime which may 
explain the elastic incommensurate magnetic response recently observed in iron-codoped Bi2201 
materials. 

PACS numbers: 71.27.+a, 73.22.Gk, 75.10.Lp 



I. INTRODUCTION 

The investigation of Hubbard type models with re- 
gard to solutions exhibiting modulated magnetic order is 
to a large extent motivated by (in)elastic neutron scat- 
tering experiments on high-temperature superconductors 
(for an overview cf. Ref. [H). Especially in the under- 
doped regime many compounds show a pronounced low 
energy spin response which is peaked away from the an- 
tiferromagnetic wave-vector. For overdoped samples of 
lanthanum cuprates it has been shown 2 that these low 
energy incommensurate spin fluctuations vanish at the 
same concentration where superconductivity disappears 
thus suggesting a close relation between both phenom- 
ena. 

There are basically three possibilities which can ac- 
count for the observed incommensurate spin response. 
First, the system can be close to a magnetic instabil- 
ity which in a random-phase approximation based treat- 
ment would strongly enhance the magnetic susceptibiliy 
at the corresponding wave-vectors. Moreover, in case of 
a superconducting (SC) d-wave gap the depletion of low 
energy spectral weight induces quasi undamped low en- 
ergy incommensurate spin excitations which merge at the 
antiferromagnetic wave- vector into the famous resonance 
peak.— ~— This scenario is very popular in optimally doped 
YBCO, since the experimentally observed energy of the 
resonance peak scales with the SC transition temperature 
T c as expected for this model. However, one could also 
envisage a situation of real symmetry breaking in the spin 
channel where (in terms of the above counting) the sec- 
ond and third possibilities correspond to ordered phases 
in the transverse and longitudinal spin channel, respec- 
tively. Since the charge couples to the square of the mag- 
netic order parameter one expects in the latter case also a 
concomitant charge modulation. In lanthanum cuprates 
(LCO), doped with Ba or codoped with Nd or Eu, « 
the third scenario is now well established. In these com- 
pounds the spin response even becomes static and con- 



comitant charge order appears, which first was evidenced 
through the coupling to the lattice^ but recently more di- 
rectly through soft resonant x-ray scatteringi&ill The fact 
that the periodicity of the charge modulation is half that 
of the spin strongly suggests that these materials have 
a so called 'striped ground state', i.e. quasi one dimen- 
sional antiphase domain walls of the antiferromagnetic 
order host the doped charge carriers. A polarized neu- 
tron scattering study on Lai^sNdo^Sro.^CuO^iL in fact 
suggests that the magnetic order is one-dimensionally 
modulated although some two-dimensional noncollinear 
structure^ may not be excluded. 

It is worth noting that stripe textures have been pre- 
dicted as stable Hartree-Fock (HF) saddle points of the 
Hubbard modeli 3 - - — before they were found experimen- 
tally in cuprates (and nickelates). However, it was real- 
ized early on that one has to include correlations beyond 
HF in order to correctly describe the experimental data 
in lanthanum cuprates. In the past years two of us have 
performed detailed investigations within the Gutzwiller 
approximation (GA) supplemented with Gaussian fluc- 
tuations which allowed for the explanation of the dop- 
ing dependent incommensurability and various transport 
propertie o 16 i , the optical conductivitj^ and magnetic 
excitation a 19 ' 20 on the basis of striped ground states. 

Non codoped lanthanum cuprates show a low energy 
inelastic incommensurability with a doping dependence 
rather similar to that of the codoped compounds sug- 
gesting the existence of some form of fluctuating and 
(or) disordered stripe order. Moreover, below a doping 
S = 0.055 the magnetic response rotates from vertical 
(i.e. along the Cu-0 bond direction) to diagonal, be- 
coming static and one-dimensional with the associated 
modulation along the orthorombic 6*-axis, also in agree- 
ment with the stripe picture. On the other hand, direct 
probes of charge order are more difficult to obtain, leav- 
ing some room for the existence of incommensurate mag- 
netic structures without charge order such as spirals. Al- 
though there are good hints from local probes like NMR, 
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NQR 21 ~— and tunneling^ 5 - that point to charge ordered 
states, it is worth examining the possibility that trans- 
verse spin (spiral) textures exist. 

Even more controversial is the situation in YBCO, the 
second class of high-T c compounds which have been in- 
tensively investigated by neutron scattering (NS) exper- 
iments.— Since in optimally doped YBCO the low en- 
ergy incommensurate magnetic response shows up below 
T C! one line of thought is that the signal corresponds to 
the dispersion of a bound exciton which is formed inside 
a d-wave superconducting gap.— ~— However, it should 
be noted that for underdoped YBCO the observation of 
an incommensurate (and even static) spin response also 
above T ,j 31 ' 32 raises questions about this interpretation. 
In any case, the lack of experimental evidence for charge 
modulation in YBCO (which if dynamic or disordered 
is of course hard to detect) makes it difficult to conclu- 
sively attribute the incommensurate magnetic response 
to stripes as in codoped LCO, though it has been suc- 
cessfully described within a model based on slowly fluctu- 
ating (or disordered) charge stripes^ 3 - Also the quantum 
oscillations found in underdoped YBCO^i are compatible 
with a reconstruction of the Fermi surface due to stripe 
orden 35 i 36 

Due to the lack of large single crystals other cuprate 
superconductors are less intensively studied by NS than 
the aforementioned compounds. However, a recent 
study has revealed elastic incommensurate magnetic 
peaks in an overdoped Bii.75Pbo.35Sri.gCu06+ 2 sample 
[(Bi,Pb)2201] codoped with iron.— Within the error bars 
the doping 5 0.23 and the measured incommensu- 
rability e « 0.21 seem to extend the relation e ~ 5 
(which holds for underdoped lanthanum cuprates^ 8 -) to 
large doping without saturation at 5 ~ 0.12. 

An alternative to stripes in order to account for the in- 
commensurate spin response in cuprates is based on the 
formation of spirals. These textures are characterized by 
a homogeneous periodic and planar rotation of the Cu 
moments with no associated charge modulation. Thus at 
first glance spirals seem to be a promising candidate for 
those cuprates in which incommensurate magnetic fluc- 
tuations have been observed without charge fluctuations, 
and corresponding theories have been recently put for- 
ward in Refs. Ei-Sl 

Within Hubbard-type models spiral solutions have 
been investigated on the basis of HFj 4 ^— slave-boson (or 
GA ) 42 ' 43 ' 45 — and dynamical mean- field 48 calculations. 
The energies of stripes and spirals have been computed 
within HF 49 > and slave-boson method s 45 ' 46 based on the 
Hubbard model. These latter investigations have found 
a strong influence of the ratio between next-nearest and 
nearest neighbor hopping t' jt on the respective stabil- 
ity of both textures. Here our starting point is simi- 
lar, however, we show that the low doping phase sep- 
aration instability of spirals can lead to the formation 
of nanoscale charge and spin inhomogeneities with sub- 



stantial spin canting. Stripes and spirals should therefore 
be viewed only as limiting cases of these more complex 
spin and charge textures. We discuss our findings with 
regard to recent neutron scattering experiments on iron- 
codoped lanthanum— and bismuth 3 ^ cuprate supercon- 
ductors. The paper is organized as follows. In Sec. II we 
introduce the formalism and present our results in Sec. 
III. Discussion and conclusions are presented in Sec. IV. 



II. FORMALISM 

Our starting point is the one-band Hubbard model 

H =J2 *« C U C J> + U Y1 n a n a> (!) 

where Ci i<T (cJ CT ) destroys (creates) an electron with spin 

a at site i, and n^ a = <\ c Ci,<y- U is the on-site Hub- 
bard repulsion and denotes the hopping parameter 
between sites i and j. We restrict to hopping between 
nearest (~ t) and next-nearest (~ t')neighbors. Our ap- 
proach is based on the Gutzwiller variational wave func- 
tion \9 g ) = P g \SD) where P g is the Gutzwiller projec- 
tor and \SD) a Slater determinant. For \SD) we use 
a state with arbitrary charge and spin order, including 
spin canting. We define the associated one-body den- 
sity as Pij ,a2 = (SD\Cj a2 Ci ai \SD). The wave-function 
optimization problem leads to a generalized Gutzwiller 
approximation 51 ' which on the saddle-point level is equiv- 
alent to the Kotliar-Ruckenstein slave-boson approach 52 . 
The derivation of the spin-rotational invariant Gutzwiller 
energy functional can be found in Ref. [53| 

E GA =^t ii (* 1 t ZiZj * j ) + A- (2) 

i,j i 

Here we have defined the spinor operators 

*, ( = (4 t ,i) 

and the z-matrix 

(Ztf cos 2 ^ + Zii sin 2 ^ [ztf - z 4 ] cos tpi 

1^- [ztx - Zi\] cos (p l Ztf sin 2 ^ + cos 2 ^ 

and for clarity spin expectation values in the Slater 
determinant are denoted by S? = p u , S i — p u , 
Sf = (plf r p^)/2, and Pll = pjf + 4"K In the limit 
of a vanishing rotation angle ip the z-matrix becomes di- 
agonal and the renormalization factors 
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- Pi + D t )(\ Pl + ^ _ A) + i pi _ _£L_ _ A) 



reduce to those of the standard GA. Spiral solutions are 
then computed by minimizing E GA with respect to a 
homogeneous rotation of spins with wave-vector Q 

Sf = S* cos(QR l ) 

Sf = S'osinCQR,). (3) 

In the present paper we usually measure modulations 
with respect to AF order Qaf = (tt, tt) (lattice constant 
a = 1) and set Q = Qaf — q- Stripe solutions are ob- 
tained by restricting the magnetization to the z direction 
resulting in a modulation of Sf with Q and a simulta- 
neous modulation of pi with 2Q similar to our previous 
work , 54 ' 55 However, we also consider stripes and other 
charge ordered states with spin canting in Sec. IIII Bl 

III. RESULTS 

In order to fix the value for the onsite repulsion U we 
refer to a previous paper where we have shown that a 
time-dependent extension of the GA with U/t = 8 can 
accurately reproduce the magnon excitations of undoped 
as revealed by neutron scattering, 56 Since U/t 
should not vary among the cuprate materials we restrict 
to U/t = 8 but investigate the dependence on the next- 
nearest neighbor hopping t'/t which from LDA compu- 
tations has been shown to specify the various high-T c 
families, 5 - Note that the results of Refs. 45.46 have been 
obtained with a significantly larger value of U/t = 12. 

A. Spirals vs. stripe states 

Fig. Q] shows the energy landscape for spiral solutions 
(hole doping measured from half-filling 5 = 0.2) for dif- 
ferent values of t'/t. With increasing \t' /t\ one observes a 
shift of the minimum from the vertical (q x , 0) and (0, q y ) 
directions towards the diagonal. Thus in the parameter 
regime relevant for LCO compounds \t'/t\ ~ 0.1... 0.2 
the incommensurate scattering direction for the corre- 
sponding doping is correctly reproduced. 

We note that also for other dopings only vertical or 
diagonal spiral textures correspond to minima of the en- 
ergy landscape in agreement with previous investigations 
i 42 1 43,45- - 47 ^^jg can easily understood from the nest- 
ing curves^ which are also shown in Fig. [T] below the 
energy contours. Basically these are q = 2kp plots (kp 
being the Fermi momentum) of the paramagnetic sys- 
tem backfolded in the < <\ x ,y < ^ quadrant. As dis- 
cussed in Ref. [58| the dominant instabilities occur for 



those wave vectors which correspond to a crossing of two 
nesting curves ('double nesting'). These double nesting 
points determine the maximum susceptibilities and are 
always found to lie along high-symmetry directions. The 
change from vertical to diagonal spirals upon increasing 
t'/t can therefore be understood from the appearance of 
an additional 'antinodal' double nesting vector along the 
(0, 0) —¥ (7r, 7r) direction which starts to dominate over 
those along the — > (n, 0), (0, tt) directions. 

The spiral is a realization of the corresponding instabil- 
ities, but the wave- vectors of energy minima and double 
nesting in Fig. [T]do not exactly coincide, as may be seen 
by comparing the nesting curves with the contour maps 
in lower parts of Fig. QJi,b,c. The Stoner criterion gives 
the correct nesting vector at the instability threshold, but 
as U/t increases the optimal q can shift. From [T] it can 
be seen that for U/t — 8 the spiral minimum coincides 
with the nesting curve near t'/t = —0.2, but is closer to 
(further from) (tt,tt) for large (smaller) t'/t. 

Fig. [5] compares the doping evolution of vertical and 
diagonal spiral wave- vectors, for t'/t = —0.2 and t'/t 

. - 1 . with the corresponding values for vertical stripes 
(cf. also Ref. [5459). The incommensurability is defined 
in terms of the wave-vector q of the dominant Fourier 
component of the magnetization as 

I Qaf - q| 

For stripes we restrict to regular, periodic stripe ground 
states with integer d (in terms of the lattice constant) 
separation between the charge stripes which amounts to 
a magnetic periodicity 2d and an incommensurability, 

1 

£ ~ 2d 

which develops a staircase structure, Fig. 2. The steps 
occur at the crossing of the corresponding energy curves 
(cf. Fig. [3]). It is also convenient to define the number 
of doped holes per unit length along the stripe v = 8d 
which is related to the incommensurability by 



Given a certain number of holes one can ask if it is 
more favorable to add them to a small number of stripes 
with a large stripe filling or in a large number of stripes 
with a small stripe filling. For widely separated stripes 
one can determine an optimum linear density v Q from 
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the minimum of the energy per hole with respect to half- 



Here E^ h is the total energy of the system doped with Nh 
holes. We find it convenient to use eh to characterize the 
competition among the different phases at all dopings. 

Fig. [3] compares eh for stripes and spirals as a function 
of doping for t'/t = —0.2 and t'/t — —0.4, respectively. 
For stripes the individual curves correspond to vertically 
oriented domain walls and are labeled by the charge pe- 
riodicity d. The optimum filling for large d is given by 
v Q = dS with S corresponding to the doping at the 
curve minimum. Notice that v a becomes independent 
of d for large d and takes the value v a = 0.55(0.4) for 
t'/t = — 0.2(— 0.4). The change of v Q with t' is expected 
according to the results of Ref. [H. 

A constant v Q implies via Eq. Q a linear relation be- 
tween doping and incommensurability. In other words 
the optimum filling determines the slope in the famous 
Yamada plot^ 

Where does the linear relation breaks down? The 
charged core of the stripe has a characteristic width £ 
so that when the charge periodicity d is larger than £ 
there are negligible interstripe interactions and doping 
proceeds by increasing the number of stripes. Therefore 
the relation is linear for d > £. In this regime the value 
of eh at the minimum is independent of doping. 

For d < £ stripe overlap becomes important and doping 
proceeds by increasing the charge of stripes. From Fig. [5] 
we see that £ ~ 4(5) for t'/t = -0.2(-0.4) which can be 
checked directly from the charge profile (see Ref. l54l ) . 

As mentioned above the incommensurability for stripes 
shows a staircase structure. For small doping (d > £), 
since interstripe interactions are negligible, one can pro- 
duce a practically continuous curve by considering com- 
binations of solutions with periodicity d and d + 1. This 
becomes very costly for d < £ which is consistent with the 
tendency of the incommensurability to develop a plateau 
at doping larger than i/ ^ ~ 1/8. 

As remarked before ) 16 i 20 ' 54 the linear relation for small 
doping and the plateau for large doping reproduce the 
doping dependence of the incommensurate low energy 
magnetic response of lanthanum cuprates as revealed 
by neutron scattering experiments^ and reproduced in 
Fig.H 

For spirals, a minimum of eh at doping 5 S indicates 
that if the frustration due to the long-range Coulomb 
interaction is ignored^, the system will lower its energy 
by phase separating in undoped regions and spiral regions 
of doping S s . Contrary to the stripes, for spirals one 
observes a continuous increase of e(S). Nevertheless, both 
phases show similar evolution of incommensurability with 
doping and indeed fairly similar values of e(S) up to S ~ 

Vol- 

For larger doping this agreement is lost due to the ten- 
dency of stripes to make wide plateaus. In fact, both spi- 



rals and stripes originate from the same magnetic insta- 
bility, the doping dependence of which has been analyzed 
for cuprate parameters in Ref. l58l . In particular, it is the 
plateau in the spin susceptibility close to Qaf = (tt, t) 
(cf. Fig. 1 in Ref. l58l ) which drives the system unstable 
towards spiral or stripe order for sufficiently large on-site 
repulsion U/t. We will come back to this point in the 
next section when discussing recent elastic neutron scat- 
tering experiments on overdoped Fe-LSCQ^S, where the 
resulting incommensurabilities are shown as solid circles 
in Fig. H 

We proceed by comparing the doping dependence of 
the energy minima which determine the orientation and 
respective stability of spirals and stripes. It should be 
noted that at low doping diagonal stripes are almost acci- 
dentally degenerate with vertical stripes with a negligibly 
preference for diagonal textures J^- Therefore other mech- 
anisms like long-range Coulomb interactions, multiband 
effects, and lattice distortions are decisive in determining 
the exact stripe state in the very underdoped regime. In 
the following we restrict for simplicity to vertical stripes. 
The envelope (dashed-dotted line) for this set of curves 
represents the doping evolution for the stripe minimum 
energy and determines the doping dependence of the in- 
commensurability as shown in Fig. [2j 

The energy minima of the spiral landscape (cf. Fig. 
[U are also shown in Fig. [3J as solid and dashed lines for 
diagonal and vertical spirals, respectively. As anticipated 
from Fig. [T]the crossover from diagonal to vertical spirals 
shifts towards higher doping with increasing \t'/t\ and 
it is found that for t'/t = —0.4 diagonal solutions are 
lower in energy over practically the whole doping range. 
However, for both diagonal and vertical spirals eh has a 
minimum at 5 S which means that for doping smaller than 
S s spiral ground states are unstable with respect to phase 
separation as mentioned above. The resulting energy of 
the phase separated solutions calculated from a Maxwell 
construction is shown as a thin horizontal solid line in 
Fig.H 

Also shown in Fig. [3] (full circle) is the energy of a sin- 
gle hole which is self-trapped in a spin-polaronic state. 
Since the corresponding charge is localized within a five- 
site plaquette with small ferromagnetic polarization, e/j 
for a many-polaron state will be independent of dop- 
ing up to the point where polaron-polaron interactions 
become noticeable. For t'/t = —0.2, both spirals and 
polarons are unstable with respect to the stripe phase, 
whereas (diagonal) spirals become the stable phase for 
larger \t'/t\, consistent with the results of Ref. [5SL 



B. Charge order states with spin canting 

The above finding that low doping spirals are unstable 
towards phase separation does not require that macro- 
scopic phase separation is the ground state. There could 
exist inhomogeneous solutions with even lower energy. 
One could envisage e.g. an elliptical spiral 6 A correspond- 
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a) t'/t=0 




b) t'/t=-0.2 




c) t'/t=-0.4 




FIG. 1: (color online) Energy of spiral solutions (per lattice 
site and in units of t) as a function of spiral momentum Q — 
(<fei <&) (in units of the lattice constant a = 1). Doping S = 
0.2. The lines in the (q x , q y , 0)-plane correspond to the nesting 
curves^, obtained for the same doping in the paramagnetic 
state. 

ing to spin structure 

Sf = So cos(a) cos(QRj) 

Sf = S , sin(a)sin(QR i ). (6) 

which connects stripe and spiral states [cf. Eq. by 
tuning the eccentricity from a = to a — tt/4. As 
discussed below there are more complex structures which 
minimize the energy. 




doping 5 

FIG. 2: (color online) Doping evolution of the spiral and stripe 
wave-vectors q for t'/t = —0.2 and t'/t = —0.4 expressed 
in terms of the incommensurability. Note that the staircase 
structure of the stripe incommensurability is only explicitly 
shown for larger doping. The solid diamonds show the dop- 
ing dependent incommensurability in LCO ('Yamada plot') 
and are reproduced from Ref. 38- The solid square corre- 
sponds to the elastic neutron scattering data on iron-codoped 
(Bi,Pb)2201 from Ref. H3- Solid circles report the incommen- 
surability of Fe-LSCO from Ref. \EH 



By performing an unconstrained minimization of the 
spin-rotational GA energy functional^- Eq. ([2]) on finite 
clusters, we find that for t'/t = —0.4 the low doping so- 
lutions are diagonal stripes with significant spin canting. 
An example is shown in Fig. 0|a) and the corresponding 
energy is indicated by the full square in Fig. EJb). This 
texture is characterized by a domain wall of the antifer- 
romagnetic order which has a fractional phase change of 
the AF order parameter Ad < tt (with 9 the angle be- 
tween the staggered magnetization and the quantization 
axis) contrary to collinear stripes which have A8 = n. 
Thus instead of building up a macroscopic phase sep- 
aration between undoped AF and doped spiral regions 
(doping S w 0.11 as can be read from Fig. [3p), the sys- 
tem prefers to separate these textures at the nanoscale, 
corresponding to the modulated structure shown in Fig. 
[4ja) . The material has clearly separated into ferromag- 
netically ordered diagonal lines which are at the center 
of 3-site wide spiral stripes with an average doping of 
5 » 0.11 and undoped AF regions located between these 
stripes. 

Research on multiferroic a 63 i 64 has shown that spin 
canting produces a force on the ligand ions via the 
Dzyaloshinskii-Moriya (DM) interaction^r— which in 
some cases can lead to a uniform dipole moment. For the 
complicate textures we find the distortion is not uniform 
but can act as an experimental signature of the canted 
phases which we describe in some detail. 

If we think of the one-band Hubbard model as embed- 
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FIG. 3: (color online) Comparison of diagonal (blue solid) and 
vertical (red dashed) spirals with the energy of stripes (thin 
grey lines) for t' jt = -0.2 (a) and t' /t = -0.4 (b). The en- 
veloping curve (dash-dotted) of the stripe energies is a guide 
to the eye. Vertical dotted lines indicate the concentration 
of zero compressibility for diagonal spirals and the thin solid 
line follows from the Maxwell construction. The full circle 
corresponds to the energy of a single spin polaron and the 
square is the energy of the spin-canted diagonal stripe which 
is shown in Fig. 0] 



ded in the cuprate lattice, then the spin canting causes 
a force on the oxygen atoms, leading to a structural dis- 
tortion of the Cu02 planes. In the presence of spin-orbit 
interactions one can add to the Hamiltonian the following 
DM and elastic term: 



HdM — ^ Au n+1 / 2 ■ fn.n+l' 

n 

h e = ^||u„+i /2 | 2 



(7) 
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where u are the oxygen displacements and the DM force 



field is given by f, 



DM — 



eij x Si x Sj. Here e„ jra+ i is 



a unit vector joining nearest neighbor Cu atoms. The 
product Si x Sj is also known as the chiral vector order 
parameter^S 

For definiteness, we assume that the spins lie in the 
CuC>2 plane, as is approximately true for the cuprates. In 
this case, only the components of Si x Sj perpendicular 
to this plane are finite, and i® M lies in the plane at 
right angles to the Cu-O-Cu bond. Once fR M becomes 
finite, He + Hum is minimized by developing a finite in- 
plane oxygen displacement perpendicular to the Cu-Cu 
bond. Due to the breaking of inversion symmetry around 
the oxygen atoms, this relaxation creates a local dipole 
moment. 



FIG. 4: a) Diagonal stripe solution with spin canting. The 
(hole) charge is accumulated on the lines with ferromagnet- 
ically aligned spins (parallel to the z-direction). b) Corre- 
sponding pattern of DM spin-currents jf M which are assigned 
to the lattice by summing the in- and outgoing currents of 
the connecting bonds, i.e. jf M = • 16 x 16 system 

with 16 holes corresponding to doping x = 1/16. U/t — 8, 
t'/t = -0.4. 



In the following we use the 'DM current' (or force flux) 

(8) 



•DM | 

Hj 



i,3 " I Si x ^j\ e i,j 



in order to characterize the ground state magnetization 



pattern 



71.72 



These currents lie in the same plane as the 



spins and for the canted diagonal stripes are shown in 
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Fig. IDJb). One observes a maximal current perpendicu- 
lar to the stripes with its direction alternating from stripe 
to stripe, so that the periodicity is the same as for the 
underlying spin structure. In contrast to the elliptic spi- 
rals, Eq. ©J which have a finite net force flux, the total 
spin force vanishes for our diagonal solution with an even 
number of stripes. The structure of the texture shown in 
Fig. 2] can be decomposed in harmonics as 



a) 



n 

Sf = J^caKQWRO 



(9) 



with q = (||, ||) and we have set Q(™) = Qaf — nq. 
The solution breaks spin rotational symmetry, so the rel- 
ative weights of the Fourier components on the x and y 
magnetization directions depend on the particular solu- 
tion or quantization axis. For the choice shown in Fig. [4] 
the amplitudes are zero in the x (y) direction for n even 
(odd). The nonzero amplitudes are given by Sf « 0.24, 
Sq ~ 0.03, and 5*1 ~ 0.1. Therefore the x-component 
of the spin structure [Fig. EJa)] induces incommensurate 
correlations atQ = (7T— ff-,1" — y|) and weaker ones at 
higher harmonics, Q= (ir — |?,7T — jg), etc. In addition, 
the y-component leads to (weaker) commensurate corre- 
lations at Q = (n, 7r) and Q = (n — n — ^r), etc.. In- 
terestingly, this implies that in neutron scattering a weak 
commensurate peak will appear, a feature that may be 
hard to distinguish experimentally from phase separation 
between commensurate and incommensurate stripes. 

For smaller values of \t'/t\, stripes do not profit from 
the transverse spin degrees of freedom and the energy 
of the spin canted solution is approximately that for the 
collinear solutions (cf. square symbol in Fig. G2[a)). Nev- 
ertheless, the finding that very different patterns are so 
close in energy suggests that stripes will be very sus- 
ceptible to quenched disorder, inducing charge-spin glass 
behavior and making it difficult to obtain clear stripes 
signatures. 
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C. Checkerboard order for larger \t' /t\ 

In general, the spin canting allows the system to make 
a compromise between ferromagnetic order (which is fa- 
vored for large |t'/t | 73 ' and AF order induced by the 
local correlations. The fact that local ferromagnetism oc- 
curs for large values of the next-nearest neighbor hopping 
was also found to play a role in the stabilization of two- 
dimensional checkerboard textures with small ferron type 
clusters for large values oft' /t = —0.5^2, These structures 
optimize the gain in kinetic energy for a collinear two- 
dimensional spin arrangement and turn out to be more 
stable than (collinear) stripes in this parameter regime. 
Based on our present findings it is thus natural to inves- 
tigate the stability of spin-canted checkerboard textures. 

Fig. [5^a) shows the minimum energy solution for 16 
holes on a 16 x 16 lattice for t' /t — —0.5. Similar to our 



FIG. 5: Spin structure (a) and DM spin currents (b) for the 
two dimensional 'checkerboard' structure with spin canting. 
Regions with larger hole density are enclosed by boxes. 16 x 
16 system with 16 holes corresponding to doping x = 1/16. 
U/t = 8, t'/t = -0.5. 



previous findings in Ref. |59j it is only for such large val- 
ues of t'/t that two-dimensional textures are stabilized 
over the (spin-canted) stripe solutions. The doped holes 
are confined to 2 x 2 plaquettes which are indicated by 
the dashed squares in Fig. [5ja) . The overall pattern con- 
sists of a 'checkerboard' array of boxes each containing 
four charge plaquettes. The planar spin structure is no 
longer incommensurate but has dominant correlations at 
the AF wave-vector together with a finite ferromagnetic 
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moment. Again, the spin structure can be characterized 
by the associated "DM spin currents" which are shown 
in Fig. EJb). Here the centers of each box appear as a 
spin current sink whereas the center and corners of the 
lattice between boxes appear as sources for the spin cur- 
rent. Note that: 1) there is no net DM spin current as 
for a spiral so there will be no net dipole moment, 2) ac- 
cording to our definition the DM spin currents are useful 
to visualize the expected lattice distortions but do not 
correspond to currents of a conserved quantity. 

We have also investigated the stability of vortex lat- 
tices which have been proposed by Finei^ but found 
them always higher in energy than the checkerboard-type 
structures. 



IV. DISCUSSION AND CONCLUSIONS 

In a recent publication,— the magnetic phase diagram 
of cuprates was computed in terms of Fermi surface nest- 
ing, which determined the critical interactions U as a 
function of doping, and the leading incommensurate Q 
vectors at threshold. The present calculations work at 
fixed U above threshold, but in general find very similar 
Q values at fixed doping. Thus, for t'/t = -0.4, Fig. 3(b) 
shows a crossover from diagonal to vertical Q as doping 
S decreases. While Ref. [58| correctly predicted this local 
minimum, it missed the fact that the global minimum is 
determined by (nanoscale) phase separation. Similarly 
for t'/t = —0.2, the Q of the vertical stripes is consis- 
tent with Fermi surface nesting 58 for 6 > 1/8, but not at 
lower doping. 

Our investigations have revealed that the nature of 
symmetry-broken states in underdoped cuprates strongly 
depends on the strength of the next-nearest neighbor 
hopping. A value of t'/t, which is appropriate for lan- 
thanum cuprates 57 stabilizes striped ground states so 
that the present analysis supports our related calcula- 
tions on charge and magnetic excitations^— as well as 
transport properties^ in these compounds. Naturally 
the wave-function of such inhomogeneous textures is de- 
termined by states far from the Fermi energy so that 
one may question the present single-band analysis. How- 
ever, we have checked by explicitly evaluating spirals and 
stripes within the three-band model that the conclusions 
of the present paper stay valid but the role of t'/t is 
played by the oxygen-oxygen hopping parameter. 

In Fig. [2] we also show data from recent neutron scat- 
tering experiments 5 ^ on Fe-LSCO. In these overdoped 
samples an elastic incommensurate spin response was 
found close to the dominant nesting vectors as extracted 
from ARPES experiments. It was thus concluded that 
the induced incommensurate response signals an inherent 
instability of the itinerant charge carriers, being different 
from the low doping stripes arising from localized Cu 
spins. Within our calculations (cf. Fig. [2]) the data are 
rather close to the incommensurability curve of stripes 
obtained for t'/t = —0.2, which is the appropriate value 



for LSCO. 

Since our calculations are based on an itinerant ap- 
proach, there is no 'dichotomy' between low doping 'lo- 
calized' spin stripes and large doping itinerant ones, but 
both appear as different limits of the same model. In 
fact, at these large dopings one can also regard the stripes 
as arising from an instability of the 'nearby' (in energy) 
paramagnet similar to the spiral case as discussed in con- 
text to Fig. [TJ In other words, increasing the doping one 
can go continuously from a stripe state at low doping 
to an itinerant spin-density wave (SDW) at high doping 
which disappears at a second order transition to become 
a uniform paramagnet. 

At fixed doping one finds an analogous SDW instabil- 
ity with increasing U. For t'/t — —0.2 and S = 0.3 the 
spin density wave instability occurs at Q = (0.677T, 7r) 
and critical interaction U c /t ~ 2. We can therefore view 
the d = 3 stripes (which have a similar periodicity) as the 
realization of this instability whereas spirals for these pa- 
rameters have higher energy. Our GA calculations yield 
a charge modulation of only 8% for the d = 3 stripes at 
5 ~ 0.25 which should be considered as an upper bound 
since our calculation neglects the effect of fluctuations. It 
may be therefore experimentally difficult to distinguish 
between a stripe and spiral state in this doping range. 

For larger values of \t'/t\ as relevant for e.g. bismuthate 
high-T c cuprates we have found that spiral textures are 
lower in energy than stripes in which canting is not al- 
lowed. If this latter restriction is relaxed, stripes and 
checkerboard states with a substantial spin canting be- 
come again more stable. This can be understood from 
the fact that spirals are unstable at low doping towards 
phase separation so that stripes and checkerboard states 
can be seen as nanoscale phase separated states. 

It is interesting to remark that the more stable stripe 
at small doping is a novel kind of stripe with a fractional 
change in the phase of the order parameter as opposed to 
the usual 7T change. These states can be experimentally 
detected because they have a finite Bragg peak weight 
both at the commensurate AF wave vector Qaf = tt,tt 
and at the incommensurate positions but can also be eas- 
ily confused with phase separation among stripes and the 
commensurate antiferromagnet phase. 

As shown in Figs. [4jb) and[5^b) the spin canting in- 
duces the flow of DM spin currents along the bonds of 
the lattice which are associated with the force flux iR M . 
Due to the breaking of inversion symmetry arising from 
the oxygen displacements, such DM spin currents may 
be associated with an electric polarization P ~ f® M 2& 

For the diagonal structure shown in Fig. 0Jb) , P would 
point along the stripe but in alternate directions from 
stripe to stripe. On the other hand, in the nearly AF 
ordered regions between the domain walls the DM spin 
currents and thus the electric polarization almost vanish. 
The spin canted stripe can therefore also be considered 
as a one-dimensional 'antifcrroelectric' with a vanishing 
net polarization^ 

In the overdoped regime and for sufficiently large \t' /t\ 
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spirals do not phase separate and under certain condi- 
tions may constitute the ground state in cuprate su- 
perconductors. In this regard, the recent elastic NS 
study 3 -^ on (Bi,Pb)2201 codoped with iron is interest- 
ing since it reveals incommensurate magnetic correla- 
tions with e k, 0.21 at a doping concentration 5 0.23. 
Our present computations shown in Fig. [5] suggest that 
the elastic neutron scattering data from Ref. |3jJ are 
compatible with vertical spirals provided that the value 
of the next-nearest hopping for Bi2201 is in the range 
—0.4 < t' /t < —0.2. LDA calculations by Pavarini and 
coworkers^ yield a value of t'/t s» —0.25 for B^S^CuOe 
where our results (cf. Fig. [3]) suggest that vertical spi- 
rals are still more favorable than diagonal ones and may 
even dominate over stripes. Naturally, our computations 
which are for perfectly modulated spirals, can only be 
qualitative since the measured coherence length of the 
magnetic modulation 37 is of the order of the iron dis- 
tance. Therefore additional effects like disorder and the 
magnetic moments arising from the Fe dopants should 
also be taken into account which is, however, beyond the 
scope of the present paper. 

In the underdoped regime it is unlikely that uniform 
spirals will survive. Although for large \t'/t\ the phase 
separated spiral has nominally lower energy than the 
stripes, the macroscopically phase separated state implies 
a charge imbalance ~ O.le per Cu which would result in 



a prohibitive Coulomb cost<£°- Even without long-range 
Coulomb interaction a nanoscale phase separated state 
has lower energy. 

Clearly the spin polarons, the checkerboard states as 
well as stripes, having charge inhomogeneities, will cou- 
ple with impurities resulting in a spin and charge spin 
glass state consistent with experiments. Since this state 
will have substantial spin-canting it becomes rather se- 
mantic to distinguish it from the disordered spiral ground 
state as proposed by Sushkov.— ~— For small \t' Jt\ we ex- 
pect that short-range diagonal stripes order dominates 
in this state as observed experimentally although with a 
correlation length of the static incommensurate response 
of the same order as the periodicity^ One should there- 
fore rather think of a disordered state with meandering 
stripe which can profit from the impurity potential. 
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